Role of Symmetry in the Transport Properties of Graphene Nanoribbons under Bias 
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The intrinsic transport properties of zigzag graphene nanoribbons (ZGNRs) are investigated using 
flrst principles calculations. It is found that although all ZGNRs have similar metallic band struc- 
ture, they show distinctly different transport behaviors under bias voltages, depending on whether 
they are mirror symmetric with respect to the midplane between two edges. Asymmetric ZGNRs 
behave as conventional conductors with linear current- voltage dependence, while symmetric ZGNRs 
exhibit unexpected very small currents with the presence of a conductance gap around the Fermi 
level. This difference is revealed to arise from different coupling between the conducting subbands 
around the Fermi level, which is dependent on the symmetry of the systems. 

PACS numbers: 73.63.-b, 73.22.-f, 71.15.Mb 
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Along with recent experimental progress in prepar- 
ing a single layer of graphite \UiM •> graphene, this two- 
dimensional (2D) electronic system has attracted exten- 
sive interest owing to its unusual band structure |5j. In 
particular, graphene can be patterned via standard litho- 
graphic techniques into new quasi-one-dimensional ma- 
terials [y, |7j, graphene nanoribbons (GNRs), which have 
many properties similar to carbon nanotubes, such as 
energy gap dependence of widths and crystallographic 
orientations |7|, |8[, and ballistic transport [9]. All 
these experimental developments and unique properties 
of graphene and GNRs provide an exciting possibility 
of new nanoelectronics based on graphene and GNRs 
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A feasible approach to realize graphene-based electron- 
ics is to construct the device junctions by connecting 
GNRs with different widths and orientations JJJ, LUJ . It 
is thus necessary to systematically study the transport 
properties of GNRs as fundamental and crucial compo- 
nents of graphene-based circuits. Until now, many ef- 
forts have been made in this direction [9|, Il2l4l4j. and 
a few intriguing phenomena were predicted, such as 
zero-conductance resonance 1 1511 . half- metallic conduc- 
tion [16[ , and valley filtering [17j • Most of these theoret- 
ical studies mainly focused on the conductance of GNR 
junctions at zero bias voltage (Vbias) rather than a finite 
Vbias- However, the results of zero bias voltage some- 
times can not represent the overall transport properties 
of systems, because the electric field induced by the bias 
can change the conductance near the Fermi level [181 ]. 
Since the electronic devices always work under a finite 
bias voltage in practice, it is essential to gain better un- 
derstanding of electronic transport properties of GNRs, 
especially upon application of a finite bias voltage. 

Thus far, in most studies of transport properties of 
GNRs, all GNRs with zigzag edges (ZGNRs) are treated 
as the same type, because they have very similar metallic 
electronic properties and transmission spectra under zero 
bias voltage [9|, Il2l - ll4j . In this work, however, we will 
demonstrate that they can be classified into two groups 



that have completely different current-bias-voltage (/- 
Vbias) characteristics (but similar band structures) with 
respect to whether there exists mirror plane a [Figs. 1(a) 
and 1(b)]. Asymmetric ZGNRs behave as conventional 
conductors with one conductance quantum under bias 
voltages, whereas symmetric ZGNRs display very small 
currents because of the appearance of the conductance 
gap depending on the bias voltage. From a two-gates 
tight-binding (TB) model, this difference in transport 
properties is revealed to originate from different sym- 
metries of ZGNRs, which is responsible for the differ- 
ent coupling properties between the subbands. More- 
over, it is shown that the current through the symmetric 
ZGNRs can be remarkably enhanced by asymmetric edge 
terminations. This feature implies that the conduction 
of ZGNRs can be controlled by changing their symme- 
tries via external modulations such as transverse electric 
fields. 



We performed quantum transport calculations of 
ZGNRs by using an ab initio code, transiesta-C [19j . 
which is based on real-space, nonequilibrium Green's 
function formalism and the density-functional theory, as 
implemented in the SIESTA approach [20[ . Self-consistent 
calculations are performed with a mixing rate of 0.05, 
and the mesh cutoff of carbon atom is chosen as 100 Ry. 
Our electronic structure calculations were performed us- 
ing the ab initio SIESTA method [20[, which is based on 
the DFT within the local density approximation. In the 
calculations, structural optimizations were first carried 
out until atomic forces converged to 0.02 eV/A. Ex- 
cept where specified, the edges of ZGNRs are termi- 
nated with hydrogen (H) atoms to remove the dangling 
bonds. Our calculations confirmed that the energy of 
the spin-polarized state (zero-temperature ground state) 
is ~20 meV per edge atom lower than that of the spin- 
unpolarized state [16|, [2l|. However, the spin-polarized 
state would become unstable with respect to the spin- 
unpolarized state at finite temperature [22[ or in the 
presence of a ballistic current through the GNRs 
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In order to simulate the experimentally detectable trans- 
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FIG. 1: (Color online), (a) An asymmetric 7-ZGNR. (b) A 
symmetric 8-ZGNR, with the mirror plane (a) represented 
by the solid blue line. The dashed rectangle denotes the unit 
cell of ZGNRs. (c) /-Vbias curves of ZGNR junctions (see the 
inset) with different widths. A bias voltage, Vbias, is applied 
to the central region with the length (I) of 8 unit cells, (d), 
(e) Conductance of 7-ZGNR and 8-ZGNR under Vbias = 0.5 
V. The conductance gap, C gap , is shown as a function of Vbias 
in the inset of (e). 



port behavior under bias voltages, we focused our study 
only on the spin-unpolarized metallic state of ZGNRs. 
To reveal the physical mechanism behind the results of 
ab initio calculations, we also employed the TB Green 
function method 14|, [l8|, [23j to analyze the transmission 
spectra with a two-gates model. 



As shown in Figs. 1(a) and 1(b), iV-ZGNRs [2J can be 
divided into two groups with respect to their symmetries: 
i.e., odd(even) N corresponds to asymmetric (symmetric) 
ZGNRs. In order to investigate the transport behaviors 
of ZGNRs under bias voltages, we carry out extensive 
first-principle calculations of I-Vbi& s curves by using a 
two-probe system [shown in the inset of Fig. 1(c)], where 
left (L) and right (R) leads are semi-infinite ZGNRs. Fig- 
ure 1(c) shows typical results for ZGNRs with different 



widths when the length I of the central region (C) is 
of 8 unit cells. It is interesting to find that symmet- 
ric and asymmetric ZGNRs exhibit completely different 
transport characteristics, although they have very sim- 
ilar metallic band structures. In detail, all asymmetric 
ZGNRs (namely, N is an odd number) have linear I- 
Vbias curves with the same slope regardless of the rib- 
bon width and length I. As shown in Fig. 1(d), this 
result arises from the presence of one conductance quan- 
tum (Go) around the Fermi level. On the other hand, all 
symmetric ZGNRs (namely, N is an even number) have 
very small currents / when |Vbi as | is smaller than a crit- 
ical bias- volt age (V c ), beyond which the currents begin 
to increase. The conductance shown in Fig. 1(e) illus- 
trates that this interesting phenomenon results from the 
presence of conductance gap (C gap ) around the Fermi 
level under a certain Vbi as - Furthermore, this conduc- 
tance gap increases with increasing | Vbias | but is slightly 
smaller than e\ Vbias | [see the inset in Fig. 1(e)], and thus 
the current remains small until | Vbias I reaches V c . 

In general, distinct transport behaviors of two classes 
of ZGNRs arise from the characteristics of their electronic 
structure. However, as shown in Figs. 2(a) and 2(b), 
asymmetric and symmetric ZGNRs exhibit very similar 
band structures, which agree well with previous studies 
[16l |25|. Around the Fermi level, there exist two par- 
tially flat bands (i.e., edge states) [25|-[27j. As wavevec- 
tor k deviates from X point, these two edge states mix 
to form bonding and antibonding states (i.e., ir and 7r* 
subbands) |l5ll28J|. In addition, there exist another two 
V-type subbands around the Fermi level [see Figs. 2(a) 
and 2(b)]. Our calculations show that the energy gap 
between two F-type subbands, 2A, is inversely propor- 
tional to the width of ZGNR, consistent with previous 

TB studies [T3,[i3- 

Clearly, different transport characteristics should come 
from something beyond the band structures. We turn to 
the wave functions of the tt and 7r* subbands, which are 
the only bands near the Fermi level and thus are responsi- 
ble for the electron transmission under low bias voltages. 
As shown in Figs. 2(c) and 2(d), the wave functions of 
symmetric and asymmetric ZGNRs do exhibit different 
characteristics. In accordance with their symmetric ge- 
ometry, the 7r (tt*) subbands of symmetric ZGNRs have 
odd (even) parities under a mirror operation. While the 
subbands of asymmetric ZGNRs have no definite parity, 
due to the absence of the mirror plane a. Then, the prob- 
lem becomes why the subband parity affects the electron 
transmission under bias voltage. 

As is well known, one important effect of the bias volt- 
age Vi,ias is that the chemical potential will change along 
the transport direction. To understand the influence of 
such change and subband parity on the transport behav- 
iors, we design an ideal two-gates model [as shown in 
Fig. 3(a)] and calculate transmission spectra using the 
TB Green function method 0, El, HI . In the system, 




FIG. 2: (Color online), (a) [(b)] Band structure around the 
Fermi level, and (c)[(d)] isosurface plots of the T-point wave- 
functions of 7i and tt* subbands for 7-ZGNR [8-ZGNR]. In 
(c) and (d), dark gray (blue) and light gray (yellow) indicate 
opposite signs of the wave function. 



local gate voltages V g + and V g - are applied to the two 
gate regions separated by a buffer region, to simulate the 
change of chemical potentials [as shown in Fig. 3(b)]. 
Figure 3c shows the calculated transmission spectra of 
such a two-gate model for 8-ZGNR and 7-ZGNR with 
Vg+=—V g -=0.25 V. Evidently, they display a symmetry- 
dependent difference, which is the same as that of ZGNRs 
under a bias obtained by the ab initio transport simula- 
tion [Figs. 1(d) and 1(e)]. Then, a clear physical pic- 
ture can be drawn from this result. As illustrated in 
Fig. 3(b), the local gates shift the Fermi level in the 
gate regions, and open an energy window (the yellow re- 
gion), within which only tt subband exists in the V^+ 
region while only 7r* subband exists in the V g - region. 
For symmetric ZGNRs, the tt and 7r* subbands have op- 
posite a parity, so they can not couple with each other to 
contribute to the transmission, and then a conductance 
gap Cgap appears [the left-hand panel of Fig. 3(c)]. For 
asymmetric ZGNRs, however, without the limitation of 
a parity, a it electron in the V g + region may hop to a 
7r* state in the V g - region and contributes to the trans- 
mission, leading to the conductance of about IGo [the 
right-hand panel of Fig. 3(c)]. 

Then, the I-Vbias curves of ZGNRs [shown in Fig. 1(c)] 
can be understood with the help of the schematic band 
diagram under Vbias [Fig. 4(a)]. Here, the solid line II 
represents the quasi- Fermi level (F q ) of the system. Lines 
I and III represent the critical energies at which F-type 
subbands appear. So only it (tt*) subband exists in the 
energy range between lines II and III (I). This feature of 
band structure is important for the unique transport be- 
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FIG. 3: (Color online), (a) Schematic configuration of a two- 
gates ZGNR model, where positive (negative) gate voltage 
V g + (V g -) is applied to the red (blue) region. In our calcula- 
tions, l g (the length of gate region) and /o are taken to be 50 
and 20 unit cells, respectively, (b) Schematic band structures 
of the left lead and two gate regions (V g + =—V g -=0.25 V). 
The dashed lines represent the Fermi level of the regions, (c) 
Transmission spectra of 8-ZGNR and 7-ZGNR in the two-gate 



model with V c 



9+~ 
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haviors of GNRs, and is different from that of armchair 
carbon nanotubes where both it and 7r* subbands (cross- 
ing at the Fermi level) span the whole energy range. The 
bias voltage makes the chemical potentials of left and 
right leads (jil and jir) be separated by |eVbias|? an d 
then only the electrons in the states between jjll and 
Hr [shaded blue area in Fig. 4(a)] can contribute to 
the current [18|- When |eVbias| < A, all these incident 
electrons into the central region are in the tt state, so 
the transmission through the system depends on whether 
tt electrons can hop to 7r* state or not. For symmet- 
ric ZGNRs, the hopping integral is zero due to <r parity 
limitation, and consequently the electrons will not flow 
through the ribbon and the current is very small. For 
asymmetric ZGNRs, however, there always exists con- 
ductance near the Fermi level [Fig. 1(d)] due to non-zero 
hopping integral between tt and 7r* electrons. On the 
other hand, when |eVbias| > A ^ eV Cl F-type subbands 
start to contribute to the conductance, and thus the cur- 
rents through symmetric ZGNRs begin to increase. 

Since the significant deviation from linear i-Vbias de- 
pendence in symmetric ZGNRs (i.e., very small currents 
under a finite bias) originates from the symmetry of the 
system, it can be expected that breaking the symmetry 
may enhance the currents. To confirm this expectation, 
for original symmetric 8-ZGNR, we use hydroxy 1 group 
(OH) to terminate one edge while using H to terminate 
the other edge, and calculate its optimized structure, en- 
ergy bands, and i-Vbias curve as shown in Figs. 4(b)- 
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FIG. 4: (Color online), (a) Schematic band diagram of the 
two-probe system under Vbias- (b) The optimized geometry 
and (c) band structure of 8-ZGNR with two edges terminated 
by H and OH respectively, (d) /-Vbias curves of the fully H 
terminated 8-ZGNR (red circles), and the H and OH termi- 
nated 8-ZGNR shown in (b) (green triangles). Vbias is applied 
to the central region with the length of 10 unit cells. 



4(d). Compared with the results of fully H-terminated 
8-ZGNR [Fig. 2(b)], the energy bands do not change 
visibly, but the current under bias voltages increases re- 
markably. This happens because the asymmetric edge 
termination breaks the original mirror symmetry, and 
consequently n and 7r* subbands no longer have defi- 
nite a parities. When a bias voltage is applied, tt and 
7r* states can couple with each other and contribute to 
the conductance around the Fermi level, which leads to 
a notably increasing current. This also implies that the 
current through ZGNRs could be tuned by applying a 
transverse electric field along the direction perpendicular 
to the ribbon axes 

In conclusion, symmetric and asymmetric ZGNRs are 
found to have completely different /-Vbias characteristics 
despite the similarity of their band structure. Asymmet- 
ric ZGNRs behave as conventional conductors with one 
conductance quantum under bias voltages, whereas sym- 
metric ZGNRs show very small currents because of the 
presence of conductance gap around the Fermi level. It is 
the different symmetry of ZGNRs that leads to different 
coupling between ir and 7r* subbands and consequently 
results in distinct transport behaviors. For symmetric 
ZGNRs, we further show that breaking the symmetry of 
electronic structure can effectively increase the current 
through ZGNRs. This might be useful in the design of 
electronic devices based on ZGNRs. 
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